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ABSTRACT 

Globular clusters are dense stellar systems that have typical ages of ~ 13 billion years, 
implying that they formed during the early epochs of galaxy formation at redshifts of 
z J> 6. Massive stars in newly formed or primordial globular clusters could have played 
an important role during the epoch of cosmological reionisation (z ;> 6) as sources of 
energetic, neutral hydrogen ionising UV photons. We investigate whether or not these 
stars could have been as important in death as sources of energetic X-ray photons as 
they were during their main sequence lives. Most massive stars are expected to form 
in binaries, and an appreciable fraction of these (as much as ~ 30%) will evolve into 
X-ray luminous (Lx ~ 10 38 erg/s) high-mass X-ray binaries (HMXBs). These sources 
would have made a contribution to the X-ray background at z ;> 6. 

Using Monte Carlo models of a globular cluster, we estimate the total X-ray 
luminosity of a population of HMXBs. We compare and contrast this with the total 
UV luminosity of the massive stars during their main sequence lives. For reasonable 
estimates, we find that the bolometric luminosity of the cluster peaks at ~ 10 erg/s 
during the first few million years, but declines to ~ 10 41 erg/s after ~ 5 million years 
as the most massive stars evolve off the main sequence. From this time onwards, the 
total bolometric luminosity is dominated by HMXBs and falls gradually to ~ 10 40 
erg/s after ~ 50 million years. Assuming a power-law spectral energy distribution for 
the HMXBs, we calculate the effective number of neutral hydrogen ionisations per 
HMXB and show that HMXBs can be as important as sources of ionising radiation as 
massive stars. Finally we discuss the implications of our results for modelling galaxy 
formation at high redshift and the prospects of using globular clusters as probes of 
reionisation. 
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1 INTRODUCTION 



C\3 Hydrogen is the most abundant element in the Universe and 
it is fundamental to galaxy formation, representing the prin- 
cipal raw material from which stars form. Approximately 
16% of t he matter content of the Universe at present is bary- 
onic (cf. ISpergel et al ] l2007t ). of which about 0.34% is cold 
gas (most of which is atomic and m olecular hydrogen; cf. 
Table 1 of iFukugita fc Peebleslliool ). From these numbers 
we conclude that the bulk of cosmic hydrogen is ionised at 
present, yet there must have been a period early in the his- 
tory of the Universe when the bulk of hydrogen was neu- 
tral. This is supported by a range of observational data that 
provides strong and compelling evidence that the Universe 
underwent an "Epoch of Reionisation" that was complete 
by z ~ 6, approximately 1 billion years after the Big Bang 
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(e.g. iBecker et all l200ll : ISpergel et afll2007l ). During this 
period the cosmic abundance of neutral hydrogen declined 
dramatically, "re-ionised" by a background of ionising UV 
and X-ray radiation whose build-up was very likely linked to 
the formation of the first g enerations of stars and galaxies 
(e.g. iBarkana fc Loebll2007l ) . 

Understanding the precise nature of the sources of this 
ionising radiation background remains an important yet 
largely unsolved problem facing modellers of galaxy forma- 
tion. It is important because reionisation is expected to have 
had a dramatic impact on galaxy formation. The presence 
of a photo-ionising background can inhibit the collapse of 
baryons onto low-mass dark m atter haloes (e.g. lEfstathioul 
ll992HThoui" fc Weinberdll996h and suppress radiative cool- 
ing a nd subsequent star for mation within dark matter haloes 
(e.g. iBenson et all [2002ah . Reionisation has been invoked 
to reconcile the apparent disparity between the observed 
abundance of satellite galaxies in the Local Group with the 
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abundance inferred from simulations of galaxy halo forma- 
tion within the favou red Cold Dark Matter framework (e.g. 
iBenson et alll2002bh . In addition, there is good reason to 
believe that reionisation was also import ant in shaping the 
faint end of the luminosity function (e.g. Benson fc MadavJ 
l2003h . the clustering of galaxies (e.g. IWvithe fc Loebll2007t) . 
the low-mass end of the HI mass function, an d perhaps the 
spati al distribution of globular clusters (e.g. iMoore et al.l 
l2006h . 

However, despite its importance for galaxy formation, 
reionisation remains a largely unsolved problem because 
very little is known about the sources of re-ionising radia- 
tion. This is because it is technically challenging to observe 
directly the Universe at z > 6, and will remain so until 
the advent of next generation instruments such as LOFAR 
(e.g. lRottgerindl2003l; IZaroubi fc Silkll2005l) and JWST (e.g. 
IWindhorst et all 120061 -1 Haimanl l2008f ). Therefore much of 
our understanding of the properties of potential sources of 
reionising radiation co mes from a combination of cosmolog- 
ic al simulations (e.g. IRicotti fc Ostrikerll2004l; Sokasian et 
al. l2004h and (semi-)analytical galaxy form ation modelling 
(e.g. lBenson et a.l.ll2006l IZaroubi et al1l2007|) . whose predic- 
tions can be tested against limits inferred from the spectral 
energ y density of the radi ation background at high redshifts 
(e.g. IDiikstra et all liooih and the electron-scattering opti- 
cal depth of the cosmic microwave background r R (e.g. Shull 
& Venkatesan lgOOSl). Not unexpected l y, UV luminous mas- 
sive stars (e.g.lWyithe fc Loebll20o"l ISokasian et~alll2004 
IWise fc Abe]||200Sl) and systems in which there is accretion 
onto X-ray lum inous intermediate mass and supermassive 
black holes (e.g. IRicotti fc Ostrikerll2004 IRicotti et aill2005l; 



IZarouMelay B ) have been su ^e d as likely so ^cTs 
of the ionising radiation background. However, there are 
many unanswered questions about the nature and origin of 
these ionising sources, and without direct observations of 
the galaxy population at z J> 6 any answers we have will be 
largely speculative. 

There are fossil records of galaxy formation in the high 
redshift Universe in our own cosmic backyard. Our Galaxy 
and others in the Local Group have large populations of 
globular clusters, the bulk of wh ich are metal poor (cf. Fig- 
ure 2 of lBrodie fc Straderll2006h and old (~ 13 billion years 
old). They appear to be relatively simple systems - their 
stellar populations appear to be generally coeval, forming 
in one or more bursts, and many aspects of their dynami- 
cal evolution have been stu died in detail by dir ect A-body 
methods (see, for example. iHurlev et alj [20071 . Praagman, 
Hurley & Power 2008). 

This apparent simplicity has led globular clusters to be 
used increasingly as probes of hi gh redshift galaxy forma tion 
(see, for example, the review of iBrodie fc Straderll2006r ) be- 
cause they may, in principle, tell us about the conditions in 
which galaxy formation proceeded at early times. For ex- 
ample, the inferred ages of old metal poor globular clusters 
implies that these systems formed at a time when the Uni- 
verse was undergoing reionisation. Reionisation is expected 
to quench star formation and therefore globular cluster for- 
mation, and so it has been suggested that present-day spatial 
distribution of metal poor globular clusters around galaxies 
can be use d to measure the r e dshift of cosmological reioni- 
sation (e.g. IMoore et al.ll2006l ; iBekki fc Yahagill2006l ). 

However, globular clusters themselves may be sources 



of ionising radiation and they co uld be potenti ally impor- 
tant for cosmological reionisation. iRicottil (|2002h (hereafter 
R02) has pointed out that massive stars in primordial glob- 
ular clusters at high redshifts are extremely luminous at 
UV wavelengths, which means that they could have been 
effective ionising sources of neutral hydrogen provided the 
radiation coul d escape freely from its source (e.g. Ricotti 
& Shull |2000| ; IBenson et al.ll2006l) . Globular clusters tend 
to reside on the outskirts of galaxies at the present epoch 
and if this was the case at early times then the fraction of 
emitted ionising UV photons that can escape without be- 
ing scattered or absorbed would have been significant. R02 
estimated that massive stars in a primordial globular clus- 
ter could emit UV photons sufficiently energetic to ionise 
atomic hydrogen {E n ^ 13.6 eV) at a rate of A 7 ~ 3 x 10 53 
s _1 over the first 4 million years; this is equivalent to a mass 
of hydrogen of between ~ 1O 9 M0 an d ~ 10 10 Mp), depend ing 
on the local recombination rate (cf. IDiikstra et aDl2004h . 

Massive stars in the local Universe tend to form with 
one or more co mpanions (e.g. lRaboudlll996l: | Mason et al.l 
1 19981 : iDelgado-Donate et al.l l2004h . and this was likely to 
be the case in primordial globular clusters. Massive bina- 
ries may evolve into high mass X-ray binaries (hereafter 
HMXBs) once the more massive star collapses to form a 
compact object, either a neutron star or a black hole. This 
compact object accretes material from its companion, either 
via a stellar wind or Roche lobe overflow (RLOF) , which re- 
sults in X-ray emission. The HMXB phase is believed to 
occur soon after the first compact object is formed and its 
duration is limited by the main sequence lifetime of the sec- 
ondary, typically ~ 10 7 years. HMXBs in our Galaxy are 
observed to emit strongly in X-rays, with typical lum inosi- 
ties of L x ~ 10 35 - 10 38 erg s" 1 (cf. iLiu et aD l2006f l . and 
there is good reason to believe that HMXBs in primordial 
globular clusters at high redshifts would h ave been as lumi- 
nous, if not more so (e.g. iDrav et al.ll2006h . 

This is very interesting because it suggests that mas- 
sive stars in globular clusters could have been important as 
sources of both ionising UV radiation and X-rays during the 
epoch of cosmological radiation. Certainly, primordial glob- 
ular clusters would have formed an abundance of UV lumi- 
nous massive stars that would have had a profound impact 
on cold star-forming gas in their immediate surroundings. 
However, could primordial globular clusters have formed suf- 
ficient numbers of HMXBs to make an interesting contribu- 
tion to the ionising X-ray background during the epoch of 
reionisation? Whereas UV photons are readily absorbed by 
neutral hydrogen, X-rays are much more penetrating (the 
photo-ionising absorption cross section of neutral hydrogen 
decreases with photon energy roughly as E~ 3 ) and can es- 
cape into the inter-galactic medium. 

In this paper we consider what fraction of massive stars 
in primordial globular clusters must evolve into HMXBs for 
these sources to make a significant contribution to the X- 
ray ionising background. We address this question using a 
Monte Carlo model of a primordial globular cluster of mass 
10 6 M Q . We assume that all massive stars form in binaries 
(see, for example, |Pray|[2006T) and explore what fraction of 
these binaries must eventually form HMXBs to be an im- 
portant source of reionising photons. 

The layout of the paper is as follows; in § 12.11 we dis- 
cuss the factors that determine the fraction of massive bina- 
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ries that survive to evolve into HMXBs (/ SU r) and, provided 
they survive, their X-ray luminosities. In § 12.21 and § 12.31 
we describe our cluster model in some detail and we com- 
pute the effective rate at which neutral hydrogen could be 
ionised by UV photons from main sequence stars and by X- 
ray photons from HMXBs. Finally, in section [3] we discuss 
the implications of our results for modelling of high redshift 
galaxy formation and the use of globular clusters as probes 
of cosmological reionisation. 



2 THE IONISING POWER OF YOUNG 
GLOBULAR CLUSTERS 

In the following subsections, we consider the factors that 
regulate the formation of high mass X-ray binaries (§ I2.1|l ; 
we present the details of our globular cluster model and 
estimate the total X-ray luminosity as a function of time 
(§ 12. 2[) ; and we determine what fraction of this luminosity 
is available to ionise neutral hydrogen (§ I2.3[) . 



2.1 The survival fraction / sur 

Not all massive binaries will evolve to become HMXBs. If a 
massive binary is to survive and become a HMXB, then the 
stars in the binary must not merge during the more mas- 
sive stars' main sequence (MS) and post-MS evolution, and 
the binary must survive the supernova of the more massive 
star. Whether or not a massive binary merges during the MS 
lifetime of the more massive star is determined principally 
by the details of stellar evolution at low metallicities and 
the initial binary separation. If the binary survives without 
merging, then the issue of whether or not the binary survives 
the first supernova depends on a number of factors, includ- 
ing the precise mass of the more massive star, the fraction 
of mass lost during the supernova and the kick velocity im- 
parted during the supernova. However, binary evolution in 
globular clusters at very low metallicities is not particularly 
well understood. Examples of complicating factors include 
the fact that lower metallicity stars are expected to retain 
more of their mass because stellar winds are inefficient at 
low metallicities; therefore these stars will be more massive 
at the end of their main sequence lives and will form black 
holes rather than neutron stars (cf. Hcgcr et al. 2003b). This 
is consistent with the expectation that the ratio of black hole 
to n eutron star remnants is m uch higher at low metallicities 
(e.g. iHeger fc Wooslevll2008h . Also, less mass is lost during 
the supernova when the remnant is a black hole rather than 
a neutron star, and black holes may suffer less violent natal 
kicks. These arguments suggest that binaries that survive 
to become HMXBs in primordial globular clusters are more 
likely to contain black holes, which suggests in turn that 
they could be more X-ray luminous than typical HMXBs 
in the local Universe, which are dominated by neutron star 
systems. 

To avoid some of the uncertainty surrounding the stel- 
lar physics in primordial globular clusters we introduce a 
survival fraction / sur which indicates the likelihood of an 
individual binary becoming an HMXB. For each primordial 
binary, we determine whether or not the binary will remain 
bound when the more massive star goes supernova by esti- 
mating the fraction of mass lost in the supernova explosion; 



if the binary remains bound, then it has a probability of / sur 
that it will evolve into a HMXB. We find that typically 70% 
of high-mass binaries become unbound following the super- 
nova of the more massive star, and so at most ~ 30% of the 
initial binary population will evolve to become HMXBs. We 
note that those that survive tend to host black holes rather 
than neutron stars. It is the fraction / sur of this remaining 
~ 30% of bound binaries which go on to become HMXBs. 
This approach allows us to investigate the potential effec- 
tiveness of HXMBs as sources of reionising radiation in a 
robust manner, even if the factors that determine the value 
of /sur are not well understood, we can speculate on what 
the most important determinants of / sur are likely to be. 



2.2 Modelling HMXBs in Young Globular 
Clusters 

We wish to assess the potential importance of the HMXB 
population in globular clusters for the X-ray background 
at high redshifts, and so we require a model of the stellar 
population in a typical globular cluster. We use a Monte 
Carlo model of a globular cluster of 10 6 stars and follow the 
evolution of its population of massive stars over its first 100 
million years, through their main sequence lives and into 
the HMXB phase. The main features of our model can be 
summarised as follows; 

• The Initial M ass Func t ion. We adopt t hree d iffer- 
ent IMFs - tho se of ISalpeteri (|l955h . iKroupal l|200lf ) and 
IChabrierl (|200ll ) - with lower and upper mass cut-offs of 
O.OIMq and IOOMq respectively. The most important pre- 
requisite for a HMXB to form is that the primary is either 
a neutron star or a black hole, which sets a low er mass limit 
of M > 8M© for the primary (cf. Figure 1 of lHeger et alj 
l2003ah ? this is the threshold for neutron star formation. The 
mean number of stars with masses in excess of 8 Mq and 20 
Mq (the threshold for black hole formation) are for Salpeter 
2621 (753), for Kroupa 11126 (3174) and for the top-heavy 
Chabrier 131935 (27024). 

• Binary Formation. We assume that all massive stars 
form in binaries. Initial binary orbita l param eters are as- 
signed following the approach of iDravl l)200rJ l - companion 
masses are drawn from a uniform distribution between 0.1 
and 100 Mq and orbital periods are distributed uniformly 
in logarithm between 1 and 10 4 days. 

• Massive Star Lifetimes. Massive stars have short 
main sequence lifetimes, the duration of which are deter- 
mined principally by their metallicit y. We obtain es t imate s 
for these lifetimes us ing t he results of Ekstrom et al.l (2008), 
ISchaerer et"ail (|l993h and lMevnet fc Maederl l|2000l) for Z=0, 
0.008 and 0.02 (i.e. solar metalicity) respectively. 

• HMXB Formation. When the more massive star in 
the binary reaches the end of its main sequence lifetime, we 
assume that it goes supernova and forms either a neutron 
star or black hole. The remnant's mass - which is determined 
by the mass of the star at the end of its main sequenc e life 
- is estimated following Figure 3 of (He ger et al . 2003a|). We 
calculate revised binary parameters (i.e. semi-major axis and 
period) although there are additional factors which compli- 
cate matters (e.g. velocity kick imparted by the supernova) . 
If the system loses more than half its mass in the supernova, 
it becomes unbound and we remove it from the list of po- 
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Figure 1. Variation of Total Bolometric Luminosity with 
Time. The light dotted curve corresponds to the contribution of 
massive stars during their main sequence lives, while the dashed 
curve corresponds to the contribution of the HMXBs. The solid 
curve indicates the net luminosity. 



tential HMXB candidates; this removes ~ 70% of binaries 
from consideration if we assume a Kroupa IMF. We then 
draw a fraction / sur of the remaining ~ 30% at random and 
consider them as HMXBs. 

• HMXB Luminosities. If a binary survives and forms 
an HMXB, we estimate its X-ray luminosity. Although bi- 
nary parameters are likely to play a role in dete rmining the 
X-ray luminosity (see, for example, Dray 2006), we take a 
simpler approach, drawing luminosities from a Weibull dis- 
tribution with a peak fixed at Lx ~ 10 38 erg/s but pre- 
venting a HMXB from accreting at greater than its Ed- 
dington limit. This sets an upper limit of approximately 
L x ^ 1.26 x 1O 38 (M/M ) erg/s. We note that this is con- 
sistent with the luminosities of compact X-ray sources in 
nearby galaxies whose X-ray bina ry populations are domi- 
nated by HMXBs (cf. Figure 1 of iGilfanov et "al]|2004h : sec 
§[3] for further discussion. 

• HMXB Lifetimes. We assume that HMXBs are active 
until the companion star evolves off the main sequence and 
goes supernova. 

Note that although we investigate the sensitivity of our 
model's results to the choice of IMF (i.e. Salpeter, Kroupa or 
Chabrier), for clarity we concentrate on models that assume 
the Kroupa IMF only. 

In Figure [T] we show how the total bolometric luminos- 
ity of our model cluster varies with time, where we have 
assumed that / sur =l, i.e. all of the ~ 30% of the initial 
massive binary population that remain bound after the pri- 
mary goes supernova. During the first few million years the 
luminosity is dominated by the contribution from massive 
stars, but this declines rapidly as these massive stars evolve 
off the main sequence and in some cases become HMXBs. 
The HMXB contribution grows rapidly once the most mas- 
sive main sequence stars go supernova (after ~ 3.4 million 



years for our assumed upper mass cut-off) and dominates 
the total luminosity from ~ 5 million years onwards. 

The total bolometric luminosity for the cluster peaks 
at ~ 10 42 erg/s during the first few million years, but it 
declines rapidly and has dropped to ~ 10 41 erg/s after ~ 5 
million years. The luminosity during this period is domi- 
nated by the contribution of the massive stars during their 
MS lives. After ~ 5 million years the decline becomes more 
gradual, falling by a factor of ~ 30 over the next ~ 80 mil- 
lion years, during which time the luminosity is dominated 
by the contribution of HMXBs. From this we may conclude 
that massive stars are indeed energetic sources of radiation 
during their MS lives and in terms of their total bolometric 
luminosity they produce as much energy as HMXBs, whose 
contribution extends over a much longer period. 



2.3 Implications for Cosmological Reionisation 

The total bolometric luminosity is an interesting number, 
but what fraction of this energy is available to ionise neutral 
hydrogen? We compute the total number of ionising photons 
emitted per second by both massive stars during their main 
sequence lives and by HMXBs and show its variation with 
time over the first 100 million years of the globular cluster's 
life in Figure [5] We model the ionis ing luminosity of ma ssive 
stars using the Starburst99 code (cf. iLeitherer et al.|[l999l) for 
the specific case of a cluster of 10 6 stars, assuming that all of 
the stars were formed in an instantaneous burst and Geneva 
tracks with high mass loss for metallicities of Z=0.001, 0.008 
and 0.02 (i.e. solar). The ionising luminosity of HMXBs is 
calculated by assuming a spectral energy distribution that 
follows a simple power-law F(E) oc E~ a between lower and 
upper energy cut-offs of S m in=0.1 eV and i5 ma x=10 6 eV 
respectively. The total energy liberated during accretion is 
then 



E tot = (A/(2-a))(R 



2-a 
max 



El 



and so given the luminosity it is straightforward to compute 
the normalisation constant A for a given HMXB. We esti- 
mate the effective number of hydrogen ionising photons to 
be 



jV 7jeff = (13.6eV)" 



F(E)dE 



(1) 



where En m is the energy above which the mean free path 
of photons becomes of order the Hubble length (i.e. atomic 
hydrogen becomes transparent to hard X-rays). We obtain 
Slim by requiring that a(En m ) — (H(z) / c) /n{z) where cr(E) 
is the ionisation cross section of neutral hydrogen, H(z) is 
the Hubble parameter at redshift z, n(z) is the mean baryon 
density and c is the speed of light. For the redshifts in ques- 
tion (z > 6) we find that i?n m > IkeV, although the precise 
value of Eiim has relatively little impact on Af T . Equation [1] 
makes the simple assumption that an energetic photon can 
ionise multiple hydrogen atoms, and so we treat secondary 
electrons that ionise hydrogen atoms as effective photons. 

In Figure [2] we show how the rate of emission of ionis- 
ing photons from massive stars (dotted curve), from HMXBs 
(dashed curve) and from ionising sources regardless of their 
nature (solid curve) varies with time, assuming F(E) oc 
E~ , which allows us to gauge the effectiveness of HMXBs 
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Figure 2. Variation of Rate of Emission of Hydrogen Ion- 
ising Photons with Time. 
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Figure 3. Ionsing Power and Dependence on Assumed 
Hardness of Energy Spectrum. 



as ionising sources as a function of f SUT . At their peak, mas- 
sive stars produce hydrogen ionising photons at a rate of 
~ 4 x 10 52 s~ 1 during the first few million years of the glob- 
ular clusters life, but this declines rapidly. We note t hat the 
peak v alue is a factor of ~ 10 lower than is quoted bv lRicottH 
(2002) but this reflects different assumptions he made about, 
for example, metallicity (Z ~ Zq) and the IMF (Salpter). 
As one expects, the effective number of ionising photons pro- 
duced by HMXBs depends strongly on what one assumes for 
/ sur and as / sur decreases, so too does the strength and du- 
ration of the HMXB population as a source of (effective) 
ionising photons. For survival fractions / sur between 50% 



and 100% the peak rate at which (effective) ionising pho- 
tons are emitted is between ~ 5 x 10 51 s _1 and 10 s _1 , but 
this declines gradually and after 30 million years the rate is 

~ w B1 s -\ 

It is important to note that we have assumed a partic- 
ular hardness for our energy spectrum (a — 1); in Figure [3] 
we show how our results depend on a, with lower values of 
a corresponding to harder energy spectra. As energy spec- 
tra become harder (i.e. as a decreases) HMXBs become less 
effective as ionising sources; this is unsurprising because the 
proportion of energetic photons increases as the hardness of 
the source increases, with a corresponding decrease in the 
number that can be absorbed by neutral hydrogen. This 
provides an interesting additional constraint on the precise 
nature of HMXBs in primordial globular clusters, if they are 
to be considered as important contributors to cosmological 
reionsiation. 

Figure [2] suggests that the hydrogen ionising power of 
a single young globular cluster is at its most effective when 
its massive star population is still on the main sequence. 
However, the ionisation cross section of neutral hydrogen is 
a strong function of photon energy, decreasing roughly as 
E~ 3 , so that the mean free path of a UV photon is much 
shorter than that of an X-ray photon. Therefore we expect 
UV photons to be most effective ionising the relatively high 
density surroundings of galaxies while X-ray photons can es- 
cape unscathed fro m the galaxy to poten tially ionise a much 
larger volume (e.g. lRipamonti et alj|2008h . Because UV pho- 
tons are absorbed in higher density surroundings, we expect 
recombination to be important and so several UV photons 
may be required t o ionise a sing l e atom of hydrogen. This 
has been noted bv lDiikstra et al. ((2004), who estimate that 
~ 10 UV photons are required to ionise a single hydrogen 
atom, compared to <J 1 X-ray photon. This is very interest- 
ing because it implies that the HMXB population is likely 
to be as important as massive stars as a source of hydrogen 
ionising radiation. Inspection of Figure [2] shows that even 
if only a few UV photo ns (rather than the ~ 10 estimated 
bv lDiikstra et all (|200 4')) are required to ionise a single hy- 
drogen atom, the relative importance of HMXBs is boosted 
dramatically and these systems now dominate as sources of 
ionising radiation. Provided the survival fraction / sur J> 0.5 
and a ~ 1, we would expect a single globular cluster of 10 6 
stars to ionise of order ~ 5 x 10 7 Mq of neutral hydrogen 
during its first ~ 100 million years. How precisely one inter- 
prets this number is not straightforward because UV pho- 
tons are absorbed locally whereas X-ray photons contribute 
to a global X-ray background. Nevertheless, it suggests that 
the HMXB population that we expect to be present in newly 
formed globular clusters at high redshifts could have been 
as i mportant as massive stars for cosmological reionisation 
fcf. lRicottill2002h . 



3 CONCLUSIONS 

Globular clusters are old, dense and relatively simple stel- 
lar systems whose dynamical properties and evolution are 
well understood. These properties have led to increasing in- 
terest in globular clusters as probes of the conditions under 
which galaxies form ed at high reds hifts (cf . iBrodie fc Straderl 
l200eh . The study of iRicottil (|2002h has suggested that mas- 
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sive stars in primordial globular clusters could have made an 
important contribution to cosmological reionisation as ex- 
tremely luminous sources of UV photons. We have explored 
the contribution that these massive stars could have made 
to cosmological reionisation once they evolved off the main 
sequence, by assuming that a fraction formed in binaries 
that evolved into HMXBs. 

Using a Monte Carlo model of a primordial globular 
cluster, we have investigated the conditions under which a 
sufficient number of HMXBs form to have a meaningful im- 
pact on the ionising power of the cluster. We assume that 
all massive stars formed in binaries and we consider only 
those binaries that remain bound once the more massive 
star goes supernova, typically ~ 30% of massive binaries. 
Of this ~ 30%, we assume that a fraction / sur survives to 
form HMXBs. Note that this is a conservative estimate be- 
cause we have neglected the effect of mass loss over the MS 
lifetime of the more massive star. 

Our results show that we require / aU r J> 0.5 if HMXBs 
are to be as effective as MS massive stars as sources of ion- 
ising radiation. Assuming / sur ~ 1 and a typical X-ray lu- 



are likely to be local analogues of primordial globular clus- 
ters, and suggests that our assumed X-ray luminosities are 
reasonable. 

It is less straightforward to assess whether or not 
HMXBs form in the numbers we might expect based on our 
model in young star clusters in "The Antennae" or "The 
Cartwheel" . The luminosity functions presented in Gilfanov 
et al. <|2004h and lWolter fe Trinchieril [|2004h can provide us 
with some insight, however. 

We find that approximately 70% of massive binaries 
are disrupted before they can become HMXBs in our model. 
Therefore we expect at most 786 HMXBs for a Salpeter IMF 
and 952 HMXBs for a Kroupa IMF to survive in a le6 solar 
mass cluster, and we introduce a factor f 3ur ^ 1 to param- 
eterise the uncertainty as to what fraction of this surviving 
population evolve into HMXBs. In the most optimistic case, 
/sur = 1. In this case, for our adopted probability distri- 
bution, approximately 10% will have luminosities in excess 
of Lx ~ 10 38 erg/s - corresponding to between 80 and 100 
HMXBs in a single cluster. These numbers will scale in pro- 
portion to fsur, so f 3Ur =0.1 implies between 8 and 10 



minosity of Lx ~ 10 erg/s, we find that HMXBs produce HMXBs more luminous than Lx ~ 10 erg/s 



ionising photons at an effective rate of > 10 s~ over the 
first 40 million years of the cluster's life. This is comparable 
to the rate at which massive stars produce ionising UV pho- 
tons during the first few million years. Note however that 
result is sensitive to the assumed hardness of the energy 
spectrum of the HMXBs - the harder the energy spectrum 
(i.e. the smaller the power-law exponent a), the less effective 
HMXBs are as ionising sources. 

By its nature, our modelling is speculative - we know 
relatively little about the conditions in primordial globular 
clusters, and little about binary evolution at low metallic- 
ities. Yet, we know that many globular clusters formed at 
high redshifts - at least as many as we can observe in the 
present day Universe - and we have a good understanding 
of the initial mass function and metallicities of the stars. 
This allows us to estimate how many massive stars there 
might have been and what their typical lifetimes were. Our 
analysis indicates that under the right conditions, HMXBs 
in primordial globular clusters can be effective sources of 
ionising radiation. If these conditions are satisfied (high sur- 
vival rate / sur , softer energy spectra) then it reveals that 
accretion power onto neutron stars and stellar mass black 
holes can make a contribution to reionisation. 

We have assumed that the typical X-ray luminosity 
of an HMXB in a primordial globular cluster is Lx ~ 
10 38 erg/s. This is more luminous than is typical of local 
HMXBs in our Galaxy and the neighbouring Large and 
Small Magellanic Clouds (e.g. the re is a single source with 
L x ~ 10 38 erg/s in the SMC; cf. IShtvkovskiv fc Gilfanovl 
2005), but these environments are atypical of the kind we 
might expect to find primordial globular clusters forming in. 
Instead, if we consider HMXB populations in galaxies that 
have undergone recent merging activity and/or show high 
star formation rate s, we find that Lx ~ 10 38 erg/s is typi- 
cal (cf. Figure 1 of lGilfanov et al.ll2004h . For example, the 
HMXB popula tions in recent galaxy mergers such as "The 
An tennae" f cf . iFabbiano et ahlhoOll ) and "The Cartwheel" 



(cf. lWolter fc Trinchierill2004h have X-ray luminosities 
cess of Lx ~ 10 38 erg/s. This is interesting because, as we 
discuss below, young star clusters in such galaxy mergers 



Of order 50 HMXB candidates are known in "The An- 
tennae". Without knowing the detailed spatial distribution 
of these candidates, it is difficult to do a straightforward 
comparison - the sources may be associated with a single 
cluster, in which case f 3ur = 1 would seem to be favoured, 
or they could be associated with 10 clusters in which case 
fsur = 0.1 would be favoured. Clearly, this is a very uncer- 
tain game to play and would require more detailed modelling 
(e.g. taking observed distribution of SSCs in the Antennae, 
populating them with HMXBs and making mock observa- 
tions), but the numbers do not seem unreasonable. 

The crucial factor here is fsur- For HMXBs in pri- 
mordial globular clusters to be an important source of 
ionising radiation, we require f sur to be of order unity. 
We expect f sur to depend on many factors, one of which 
will be metallicity; f sur may decrease as stars become 
more metal rich, and so it may very well be the case 
that fsur = 0.1 rather than ~ 1 in a system such as 
the Antennae. Establishing how different factors impact 
upon fsur and how these factors may depend on redshift 
is something that we would like to investigate in future 
work, but at present the discussion is necessarily speculative. 

We have focussed on the ionising power of HMXBs in 
primordial globular clusters, but it is important to note that 
the longer mean free path of X-rays compared to UV means 
that HMXBs in galactic discs could have been as effective as 
HMXBs in globular clusters. We focus on globular clusters 
because we have a reasonable understanding of what the 
IMF could have been, and therefore we have a reasonable 
idea of how many HMXBs were likely to form. In contrast, 
we know very little about the initial mass function in star- 
forming galactic discs at high redshifts, but this is not to 
say that HMXBs in discs should be discounted. Rather it 
suggests that stellar mass black holes could make an inter- 
esting contribution to reionisation and merit further study. 
The contribution of globular cluster HMXBs to the reionis- 
ing background considered in this paper could be considered 
as a lower limit. 

More broadly our analysis strongly indicates that glob- 
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ular clusters could have been important ionising sources in 
the high redshift Universe. This has interesting implications 
for galaxy formation at high redshifts and the use of globu- 
lar clusters as probes of reionisation. As yet, we do not have 
a satisfactory theory for the formation of globular clusters 
within a cosmological context. Nevertheless, if they were to 
be effective sources of UV radiation then they must have 
moved quickly from their formation site, presumably a co- 
coon rich in cold dense gas, to the rarified environs of the hot 
gaseous halo; if this did not happen on a timescale shorter 
than a fraction of a few million years, a good fraction of their 
UV radiation would have been absorbed by the cold dense 
gas. This would favour globular cluster forming in gas-rich 
mergers akin to the "super star clusters" that we observe 
in present da y mergers (e.g. IWhitmore fe Schweizerl 1 19951 : 
lMeurerlll995l) rather than in the gas-rich discs fe.g Kravtsov 
& Gnedin l2005ir ~ 

We have argued that ionising radiation from the glob- 
ular clusters could have been effective in ionising neutral 
hydrogen surrounding their host galaxy, but this radiation 
could also damage cold dense gas within the galaxy (both 
neutral and molecular), thus affecting the star formation 
rate and potentially suppressing further star formation. This 
has particular implications for the use of globular clusters 
as probes of cosmological reionisation. The argument is that 
the spatial distribution of globular clusters around galaxies 
and in galaxy clusters can be used to measure the epoch of 
reionisation - the more centrally concentrated the distribu- 
tion of globular clusters, the earlier the epoch of reionisation. 
Yet if the star formation efficiency in a galaxy is high, and 
consequently the rate at which globular clusters form is high, 
then we might expect feedback from early-forming globular 
clusters to be ex tremely damaging f or later-forming globu- 
lar clusters (e.g. iMoore et al.ll2006h . This would suppress 
the star formation and globular cluster formation efficiency. 
How could one separate this local effect from the effect of a 
globally driven period of reionisation? This is a question we 
shall be pursuing in future work. 
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